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We investigate the amplification of tfie field induced by white noise. In the present 
study, we study a stochastic equation which has two parameters, the energy uj{k) of a free 
particle and the coupling strength D between the field and white noise, where the quantity 
k represents the momentum of a free particle. This equation is reduced to the equation with 

one parameter a(k) which is defined as a{k) = D \i-o[k) j . We obtain the expression of 

the exponent statistically averaged over the unit time and derive an approximate expression 
of it. In addition, the exponent is obtained numerically by solving the stochastic equation. 
We find that the amplification increases with a{k). This indicates that white noise can 
amplify the fields for soft modes if the mass m of the field is sufficiently light and if the 
strength of the coupling between the field and white noise is sufficiently strong, when the 

energy u{k) is equal to \/ m? + k'^. We show that the a{k) dependence of the exponent 
statistically averaged is qualitatively similar to that of the exponent obtained by solving the 
stochastic equation numerically, and that these two exponents for the small value of a(k) 
axe quantitatively similar. 

§1. Introduction 

Over the past few decades, a considerable number of studies have been made 
on the phenomena induced by noise and it is understood that noise plays impor- 
tant roles in many branches of physics. Well-known phenomena are stochastic reso- 
nance^)"^^ , stochastic synchronization,^)'^) noise induced propagation,^) phase tran- 
sition induced by multiplicative noise,®) noise enhanced phase locking,^) coherence 
resonance,^'') and so on. These examples indicate that noise plays essential roles in 
some systems. Therefore, the inclusion of noise is important to explain some physical 
phenomena. A phenomenon that noise may be important is the particle production 
that thermalize the system as discussed in the study of the early universe and heavy 
ion collisions. The number of produced particles may be enhanced by noise, because 
noise can supply energy. Therefore, it is worth to investigate the amplification of 
the field by noise. 

The amplification of the field has been investigated in terms of particle produc- 
tion, because the produced particles affect the time evolution of the system. An 
mechanism of the amplification is spinodal decomposition^^) that the field is am- 
plified through the roll-down from the top of the potential (false vacuum) to the 
minima of the potential (true vacuum), because the mass squared is negative around 
the top of the potential hill. Another mechanism is parametric resonance^'^)"^^) that 
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an oscillating field amplifies other fields. The effects of noise on parametric reso- 
nance have been studied. ^^^'^^^ Noise modifies the amplification quantitatively. For 
example, noise shifts enhanced modes. ^^-^ 

These mechanisms are well-known as the mechanism of the amplification of the 
field through phase transition. The amplification by spinodal decomposition occurs 
in the beginning of the phase transition and that by parametric resonance occurs in 
the end of the phase transition. Particle does not produced by parametric resonance 
if no oscillating field exists, because no parametric resonance occurs in the absence 
of an oscillating field. However, another mechanism of the amplification may work 
even in such a case if noise exists. Phenomenological equations for such processes 
are similar to the equation of a pendulum. 

The oscillator with randomly varying mass and/or friction coefficient has been 
studied in order to comprehend the effects of noise. In Ref. 26), the equations 
of motion were described with amplitude-angle variables and the condition for the 
growth of amplitude was acquired. In the similar manner, the phase transition was 
investigated for a pendulum with a randomly vibrating suspension axis.^''^ Recently, 
K. Mallick and P. Marcq investigated the nonlinear oscillator with white or colored 
noise. ^^^"^'^^ In their studies, the equations of motion were described with energy- 
angle (or action-angle) variables. The power exponents of time development were 
given by solving the (effective) Fokker-Planck equation which describes the evolution 
of slow variable. In these studies, the Fokker-Planck equation is used as a main tool. 

In this paper, we investigate the effects of noise on the amplification of the fields 
by adopting an equation with white noise. We attempt to derive an approximate 
expression of the exponent and to obtain the exponents for various values of the 
parameter. The role of physical parameters for the amplification is obviously shown 
by the evaluation of the exponents when white noise exists. In ^ a basic equation in 
the present study is introduced. The equation describes a harmonic oscillator with 
varying mass. An approximate expression of the exponent is derived by averaging 
the term with white noise statistically and by using the steepest descent method, 
without using the Fokker-Planck equation, though the Fokker-Planck equation is a 
powerful tool. In ^ we solve the stochastic equation introduced in ^numerically. 
Many trajectories can be obtained as numerical solutions and the exponents are 
extracted by averaging these trajectories. We compare the values of the exponents 
obtained from the expression derived in ^with those of the exponents obtained by 
solving the stochastic equation. Section 0] is devoted to conclusions. 

§2. Amplification induced by white noise in a scalar field theory 

2.1. Equation of Motion 

In the beginning of this section, we introduce the parameters and the variables 
used in the following discussion. The quantity Lo{k) represents the energy with 
momentum k. n{t) is a white noise which has the properties 

{n{t)) = 0, (2-la) 
{n{t)n{t')) = 5{t-t'), (2-lb) 
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where (• • • ) represents the statistical average. The starting point of this study is the 
following equation which describes the motion of the field cp for finite modes near 
the bottom of the potential with white noise: 

^ + (u;{k)) ' + Dn{t)^ = 0, (2-2) 

where the parameter D is the coupling strength between the noise and the field. 
Without loss of generality, we can set this parameter positive or zero. No oscillating 
term appears in Eq. (|2-2p . while an oscillating term appears in the coefficient of 
(p in Ref. 14), 15). Equation (|2-2|) can describe approximately the system with an 
oscillating field with quite small amplitude or without an oscillating field. We 
regard Eq. (\2-2^ as a basic equation in this study. It is not the aim to derive Eq. H2-2() 
in this paper. Following discussion is reasonable, once Eq. (|2-2p is introduced. 

In order to represent Eq. ()2-2() by dimensionless variables, we rewrite Eq. H2-2|) 
in terms of the new variable z given by z = uj{k) t, we obtain 



dz^ 



+ (l + a(A?) r(z;fc))0 = O, (2-3) 



-3/2 ^ / ^ \ 

where a{k) = Dluj{k)\ and r{z;k) = la;(fc)) n{z/uj{k)) which has the 



property 

{r{z;k)r{z';k)) = 5{z - z'). (24) 

The momentum k does not play an essential role in the calculation of the ampli- 
fication from Eqs. (|2-3j) and (|2-4|) except the momentum dependence of the field 
amplification. The amplification is determined by only the value of a(k). Equa- 
tion (|2-3|) is rewritten by introducing = d(j)/dz and dW = r{z; k)dz , because it 
is not easy to handle Eq. (|2-3|) numerically. It should be noted that dW is a wiener 
process. Therefore, Eq. (|2-3() is rewritten as 

dcj) = Pcf, dz, (2-5a) 
dp^ = -(f)dz- a{k) (p o dW. (2-5b) 

We regard Eqs. ()2-5a() and H2-5b|) as stratonovich equations, (o represents a stratonovich 
equation.) Stratonovich equation is easily converted to Ito equation. The Ito equa- 
tions corresponding to Eqs. (|2-5a() and (|2-5b|) take the same form, that is, 

d(j) = Ptf, dz, (2-6a) 
dp^ = -(j)dz- a{k) (p dW. (2-6b) 

We solve Eqs. ()2-6a|) and 1)2 -Gbl) numerically in the next section. 
2.2. Estimation of the Exponents 

Equation (|2-3|) describes the motion of a harmonic oscillator when a{k) is zero. 
Then the solution is a sine function for a{k) = 0. Hereafter, we omit the argument k 
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of a{k) when no confusion occurs. The coefficient + a{k) r{z; k)j can be negative 

for a 7^ 0, because r(z; k) is a random variable of zero mean value, {a is not negative, 
because D is larger than or equal to zero.) In this subsection, we estimate the 
magnitude of the exponent Ag which describes the increase of the amplitude A{z) 
of the oscillator in time width Az: A{z + Az) = exp{Ag)A{z). The amplitude does 
not increase if the coefficient of (p in Eq- is positive. Therefore we pay attention 
only to this coefficient which is negative. Then the exponent Ag in an extremely 
short time step Az, in which the random variable r(z; k) is constant, is given by 

Ag = AzO(^-l- a{k) r{z; k)^ (^-1 - a{k) r{z; fc)) , (2-7) 

where 0{x) is a step function which is 1 for x > and for x < 0. If the time step 
Az is large enough compared with the time length in which the random variable 
r{z; k) is constant, the statistical average of Ag should be taken. In this sense, the 
width Az is the resolution of the observation. To make this fact clear, we divide the 
time step Az into regions, in which the random noise r(z; k) is constant. When 
the amplification of the field is given by 4'{zi + {Az)/N) = exp[[3i{Az) /N)4>{zi^i) 
for the region i (The regions are distinguished by the index i.), the amplification of 
the field in Az is given by 

0(z + Az) = JJexp (^A <A(z) = exp [^^Y.^"^^)^ (^-S) 

The quantity f3i{Az) is the exponent in the case that the random noise is constant in 
time step Az. Eq. (|2-8() implies that the exponent in Az is the statistical average of 
the realized exponents Pi{Az), when the noise r(z; k) varies frequently in Az. If Pj is 
pure imaginary, it does not affect the amplification. Then step function is included 
in Eq. ()2-7|) as stated above. Therefore we take the statistical average, because the 
Gaussian white noise varies for very small time intervals. As in the previous subsec- 
tion, we define AW = r{z; k)Az. The probability distribution function P{AW) of 
the wiener process AW is a Gaussian distribution of zero mean value and of variance 
Az, that is, 

P{AW) = , ^ exp f-i^^^ . (2-9) 
We estimate the amplification by exp [{Ag)st], where iAg)st is defined by 

/oo 
d{AW)P{AW)Ag. (2-10) 
-oo 

The exponent {Ag)st with Eq. (|2-9p can be rewritten as 
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Equation (|2-1H) is easily evaluated for a ^ 1 and a ^ 1. The exponent {Ag)st for 
a S> 1 is given by 



St 



^(3/4) ,1/2,^^,3/4 
23/4^^1/2 



n3/4 



The exponent {Ag)st for a ^ 1 is given by 

/ Az\ 



St 



a 
2372 



exp 



Az 



(2-12) 



(2-13) 



We introduce an arbitrary integer n which is larger than or equal to 1 in order to 
evaluate Eq. ()2-ll|) approximately for other values of a. Rewriting Eq. ()2-lH) in 
terms of the new variable x given by t^/^ = x'^ and using the method of steepest 
descent, we finally obtain 



-1/2 



exp 



Az 
'2a 



— - f { x^"^ 



fn (x) 



in) 



X erfc 
1 



>-l/2^W 



J n 



(2) (J^) 



1/2 



^4n ^ ^A^x^n 



a 

{AzYl^ 



+ (1 — 3n) Inx, 



+ 



a 



{Azy/'^ 



a 



+ 6 



1 

2n 



(2- 14a) 
(2-14b) 

(2-14c) 



(2) 

where /A (x) is the second order derivative of fn{x) with respect to x and the error 
function is defined as 



2 [°° 

erfc (x) = --j= J dy exp (— y^) • 



(2-15) 



For example, {Ag)st from Eq. (|2-14a|) for Az ^ 1 is approximately evaluated as 
{Ag)st ~ 0.428ai/2(-^^)3/4 ^Yien n is set to 1 and {Ag)st ~ 0.400a^/2(^^)3/4 ^-^^^ 

n is set to 2. The quantity {Ag)st can be evaluated directly from Eq. (|2-11|) for 
Az < 1: {Ag)st ~ OAlla^/'^iAzf/'^. Then, it is expected that Eq. (|2-14a|) gives the 
approximate value of Eq. (|2-11() . Our purpose is to evaluate the exponent in the unit 
time Az = 1. We cannot take Az <^ 1, because the exponent (Ag)st is statistically 
calculated. Then, we put Az = 1 in the following discussions. In the next section, 
we compare Eq. ()2-14a() with the exponent obtained from the numerical solutions 
for the equation of motion. 



§3. Numerical results 



In this section, we solve Eqs. ()2-6aj) and 1)2 -Sbl) numerically for various values 
of a. The initial conditions are </>(0) = 1 and d(l){z) / dz\ = in the present 
calculations. We apply the Euler-Maruyama scheme to the stochastic equations. 
The time step in z is set to 0.05 in these numerical calculations. 
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Fig. 1. Typical time evolution of the field for 
a = 0.6. 
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Fig. 2. The time evolution of the mean value 
of for a = 0.6. One batch has 500 sam- 
ples and the mean value over 20 batches, 
<j){z), are calculated numerically . Thick 
curve represents <j){z) and vertical bars indi- 
cate the confidential interval of 50%. These 
bars are depicted every 20 points. 



Figure n displays the time evolution of (/){z) for a = 0.6. Apparently, the field 
is amplified because of white noise. One trajectory can be calculated when the 
sequence of noise is given. Thus this amplification may not occur frequently. Then 
we should calculate many trajectories and take the average of them. Therefore, we 
calculate the mean value of the trajectories of the field (j)[''\z), where the subscript 
i distinguishes different batches and the superscript (j) indicates a trajectory in a 
certain batch. One batch has 500 samples (trajectories) and the mean value Mi{z) 
of these trajectories is calculated. We take 20 batches and calculate the mean value 
of Mi{z). Then, the mean value over 20 batches, ^(-z), is given by 

1 ^° 

i=l 

500 

The averaged behavior (t){z) for a = 0.6 from z = 100 to 150 is displayed in Fig. |21 
Thick curve represents the mean value (t){z) and vertical bars indicate the confidential 
interval of 50%. These bars are depicted every 20 points in this figure. These bars 
become wide with time z, because the present process is a wiener process. It is found 
from this figure that the field is amplified by noise. We use (piz) for various values 
of a in the subsequent evaluation of the exponents. 

Figure 131 displays the exponents in the unit time dg\Az=i and the approximated 
exponents in the unit time {^g)st\Az=i ^^^^ various values of a. The exponent dg\j^z=i 
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Fig. 3. The exponents in the unit time dffl^^^j and {Ag)st\^^^j^, in (a) the narrow region (a < 3) 
and (b) the wide region {a < 20). The symbol "x" represents the data points for dg\^^^^ 
obtained from (j>[z) and the thick curve is a guide. The dotted curve represents {Ag)at\^^^-i 
obtained from Eq. (I2-11II and the dashed curve represents {Ag)Bt\^^^-^^ obtained from Eq. 12-14al l 
with n = 2. 



is estimated from <t){z) in the region oi z > 100 in order to decrease the effects of 
the initial conditions. This estimation is performed as fohows: 1) Making the set 
of (zj, In (^(zj)), where Zj is the time at which (t){zi) is local maximum and positive. 
2) Fitting the above set with a linear function. The coefficient of the time z is 
adopted as dg\Az=i- The exponent {Ag)st\Az=\ calculated in two ways. One 
[{Ag)st\Az=i (f^ll)] is obtained by using Eq. (|2-11|) and the other [[Ag)st\Az=i 
(app. : n = 2)] is by using Eq. (|2-14a() with n = 2. In Fig. |31 (a) and Fig. 0] 
(b) , the symbol "x" designates points ''^dg\/\z=i^ and the thick curve is a guide. 
The dotted curve represents "(Z\(7)st|^^=]^ (full)" and the dashed curve represents 
''^ {Ag)st\Az=i (app. : n = 2)". In both figures, the dashed and dotted curves are 
quantitatively similar. The exponent dg\Az=i is close to the exponent {Ag)iit\Az=i 
for small values of a as in Fig. |21(a) and the difference in magnitude between dg\Az=i 
and (^5)8*1^2=1 increases with a as in Fig. |21(b). Therefore, the exponent dg\Az=i 
can be approximately evaluated with the exponent {Ag)st\Az=i ™ region of small 
values of a. Obviously, numerical results imply that dg\Az=i is a monotone increasing 
function of a. In these figures, the a dependence of (^i?)st|^2=i ™^ th&t of dg\Az=i 
are similar. Note that it is inferred that the exponent dg\^^^^ estimated from (t){z) 
is equal to or larger than the exponent {Ag)st\Az=i-> because the arithmetic mean of 
positive real numbers is equal to or larger than the geometric mean of them. When 
a field 4>j{z) is given by (l)j{z) = exp(7j)(/)(0) (The index j distinguishes trajectories.) 
and the number of the trajectories is M, the following inequality is satisfied: 



10(^)1 



M ( M \ 

-5]exp(7,)|<A(0)| >exp h^E^i ^^'^^ 



8 



M. Ishihara 



The exponent of the right-hand side of Eq. (|.S-2|) is the average of exponents obtained 
from the trajectories. Assuming that this exponent is equal to the average of the 
exponents in time (Eq. H2-1U() ). we can conclude that the exponent dg\^^^^ is larger 
than or equal to the exponent (^5)8* | ^2=1- 

§4. Conclusions 

We investigate the amplification of the field induced by white noise. We obtain 
the expression of the exponent statistically averaged and the expression of the expo- 
nent approximated with the method of steepest descent. In addition, the exponents 
are extracted from the numerical solutions of the stochastic differential equations. 

We summarize the results: 1) White noise amplifies the fields, especially in the 
region of large values of a{k). This fact indicates that the amplification for soft 
modes is strong compared with that for hard modes, when the energy uj{k) is equal 

to \/ -|- k"^, where m is the mass of the field (f>. The amplification weakens with the 
mass of the amplified field. 2) The expression of the exponent statistically averaged 
is reliable to estimate the magnitude of the exponent for the small values of a{k). In 
particular, the expression approximated with the method of steepest descent is easy 
to use, because this expression is written with only well-known functions. These 
results imply that the soft modes can grow on the vacuum that is located at the 
bottom of the potential, if the coupling between the field and white noise (i.e. the 
parameter D) is sufficiently strong and if the mass of the field (j) is sufficiently light. 

As pointed in H2.21 the coefficient of (p in Eq. H2-3|) can be negative because of 
noise. The amplification in the present study is caused by the negative coefficient. 
In comparison, the amplification in spinodal decomposition process is induced by 
negative mass squared. Therefore, soft modes are amplified strongly in both the 
present process and spinodal decomposition process. 

It is expected that the stochastic equation used in the present study appears in 
various fields of physics, for example, in the study of the early universe and chiral 
symmetry rebreaking. Though the parameter a{k) may be one or less, the field 
can be amplified if a{k) is larger than 0.3. (See Fig. Ol^a).) The parameter a{k) 

decreases with the mass m if the energy uj{k) is equal to \/ m? + k"^ . Therefore, the 
amplification is strong for shallow potentials. 

In the present study, we treat the linear equation with white noise. However, 
nonlinear terms and colored noise appear in general. We would like to handle the 
effects of these terms in the future studies. 
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